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Abstract 

We solve the functional equations 



1 1 1 

f{x) f(y) f(z) 
/'(*) f\y) f(z) 



0.. 



1 1 1 

f(x) g{y) h{z) 
f{x) g'(y) h'(z) 



for suitable functions /, g and h subject to x + y + z = 0. These equations 
essentially characterise the Weierstrass p-function and its degenerations. 
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1 Introduction 

The purpose of the following note is to present a simple and direct proof of 

Theorem 1 Let f be a three-times differentiable function satisfying the func- 
tional equation 



111 

f(x) f(y) f(z) 
f\x) f'(y) f'(z) 



x + y + z = 0. 



Up to the manifest invariance 



f{x)^af{6x) + P, 

the solutions of (fl|j are one of f(x) = p(x + d), f(x) = e x or f(x) = x. Here p 
is the Weierstrass p-function and 3d is a lattice point of the p-function. 

In fact our approach gives a simple proof of 

Theorem 2 Let f , g and h be three-times differentiable functions satisfying the 
functional equation 



1 1 1 

f(x) g(y) h(z) 
f'(x) g'(y) h'(z) 



0. 



x + y + z = 0. 



(2) 



Up to the manifest invariance 

f(x) -> a/(5x + 7i) +0, g(x) -> ag(5x + 72) + 0, h(x) -> ah(5x + 73) + p, 



where 71 + 72 + 73 = ; the nonconstant solutions of (fjj are given by f(x) = 
g(x) = h(x) = e x , x, or p(x). If (say) h(z) is a constant then either 

1. One of the functions f{x) or g{y) is the same constant as h(z), in which 
case the remaining function is arbitrary, or 

2. f(x) = g(x) = e x . 

Remarks: (i) In fact the exponential and linear function solutions satisfy (|1|) 
and (H) without the constraint x + y + z = 0. 

(ii) The theorems immediately give the general analytic solutions to the same 
functional equations viewed as functions of a complex variable, showing that the 
solutions are in fact meromorphic. 
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(Hi) The arbitrary constant 5 in the invariance of ([I]) is accommodated in the 
Weierstrass jp-function solution by the homogeneity relation p(tx; t~ 4 g 2 , t~ e gs) = 
t~ 2 p(x;g 2 ,g 3 )- 

(iv) Weierstrass has shown || that any meromorphic function possessing an 
algebraic addition formula is either an elliptic function or is of the form R(z) or 
R(e Xz ), where R is a rational function. A priori the functional equation ([]]) is 
distinct from assuming / possesses an algebraic addition formula. 

Interest in (ED) arises from a question of mathematical physics: What one- 
dimensional quantum mechanical models with pair-wise interactions have ground 
states of product type? In addressing this question || the functional equation 

F(x)F(y) + F(y)F(z) + F(z)F(x) = G(x) + G(y) + G(z), (3) 

appears (with x + y + z = 0), where on physical grounds F(x) is taken to be odd. 
By applying (d x — d y ) d x d y to this equation we obtain ([!]) with f(x) = F'(x) and 
similarly (|1|) may be integrated to yield (|3]). The models that arise in the solution 
of this question include the Calogero-Moser-Sutherland models. They are rich 
in interesting mathematics involving representation theory, harmonic analysis 
and special functions (see for example ||); the classical analogues of the models 
yield completely integrable Hamiltonian systems. Now the solutions that yield 
these models were obtained assuming the function F(x), or equivalently f(x), 
to be meromorphic with at least one pole. With such an assumption it is easy 
to show that f(x) = p(x + d) is the general solution, and further requiring f(x) 
to be even dictates d = 0. Although one can in fact show that there are no 
(nonconstant) even entire solutions to ([[]), so answering the physical question, 
the general solution to ([!]) appeared difficult to obtain. Indeed, the functional 
equation (|2|) was introduced in O as a means to understand (|T|),Q but the proof 
of the main theorem of || fails to result in a direct proof of Theorem 1. Here 
we present a simple and direct proof of ([I]) that also allows a new and simpler 
proof of (0). 

2 Proofs 

The strategy of our proof is first to isolate a necessary condition for nonconstant 
solutions of (|l|) and (fj) in the form of a differential equation. The solutions 
of this differential equation are given in terms of the Weierstrass p-function 
or one of its degenerations, and the second step is to determine the various 
free parameters that arise in the solution. This will prove the theorems for 
nonconstant functions /, g and h. Finally we consider the case where some of 
these functions are constant. 

1 We are grateful to V.M. Buchstaber for informing us on this matter. 
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While the approach to solving a functional equation via an associated dif- 
ferential equation(s) is standard the simplicity of our proof depends on 
one rather nonobvious step that we wish to highlight in advance. First, by 
taking (0) and various of its derivatives we obtain an equation of the form 
F(f(x),g(y)) = involving f(x), g(y) and their derivatives. Viewing x (say 
x = xq) as fixed we have in general a nonlinear ordinary differential equa- 
tion for g(y) which may in principle be solved. While such a solution satisfies 
F(f(x ), g(y)) = 0, it needn't in general satisfy F(f(x),g(y)) = 0, the various 
derivatives d£ F(f(x), g(y))\ xo = (supposing they exist) yielding further dif- 
ferential equations for g(y). These further equations give us restrictions on the 
allowed functions g(y) and between them one may eliminate various derivatives of 
g appearing. For example, given F(f(x Q ), g(y)) = and d x F(f(x),g(y))\ XQ = 
one could could choose to eliminate the highest derivative of g appearing; sim- 
ilarly one can use further partial derivatives to eliminate additional derivatives 
of g. If we suppose that F(f(x),g(y)) = determines g(y) then the equation 
F(f(xo), g(y)) = 0, together with an appropriate number of further derivatives, 
also determines g(y)- The nonobvious step in our proof is to provide a func- 
tional F(f(x),g(y)) = that alone readily gives g(y). We will remark in the 
course of the proof when this is done, and note the several perhaps surprising 
simplifications that follow. 

Lemma 1 Let f , g and h be three-times differentiable, nonconstant functions 
that satisfy ^). Then each satisfies the (same) differential equation 

w'(x) 2 = p 3 w(x) 3 + p 2 w(x) 2 + pi w(x) + Pq. (4) 

Proof: We begin by deriving several algebraic consequences of assuming that 
the (nonconstant) functions /, g and h of (Q) are iV-times differentiable. Let 
1 < k, I, s < N. The algebraic identities we obtain yield a large supply of 
functional equations involving only the functions f(x), g(y) and their derivatives. 
We will obtain (|]) by eliminating an appropriate derivative. A minimum choice 
of N = 3 will arise in the proof. 
Set d = d y — d x and let 

a k = d k - 1 (g(y)-f(x)), b k = d k (g(y) + f(x)), c k = d k (g(y) f(x)). (5) 
Then differentiation of (0) yields N equations, 

a k h'(z) - b k h(z) + c k = 0. 
Comparing any two of these equations shows 

{a k bi- aib k )h(z) = a k ci- aic k , (a k bi - aib k )h'{z) = b k q, - bic k , (6) 
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while comparison of any three yields 



0>k 0-1 0-t 

h k b s 
Ck ci c s 



(7) 



Consider z 
equations (I 
results in 



x y as a function of x and y. Differentiating the first of 
with respect to y say, and comparing with the second equation 



a k ai a k dydi — ai d y a k 

bk k a k dybi - a t d y b k 

Ck ci a k dyCi - ai d y c k + b k ci- b { c k 

Similar expressions result upon differentiation with respect x. 

Observe that at this stage (|7|) and (||) and their linear combinations provide us 
with many functional equations involving only the functions f{x), g(y) and their 
derivatives. In particular we may eliminate various combinations of derivatives 
between them. For example, by taking k — 1, I — 2 and s = 3 the linear 
combination (||) - axx(^), 



ai a 2 atg"(y) - a 2 g'(y) - aia 3 

b\ b 2 ai g"'(y) - a 2 g"{y) - ai b 3 

ci c 2 aid y c 2 - a 2 d y ci - (X1C3 + (b\c 2 -b 2 ci) 



0. 



(9) 



cancels each of the g"'(y) derivatives appearing in the third column. This ex- 
pression, quadratic in g"(y), factorises to give 







/'(*) g\y) - g'{yf - f(x) g"{y) + g(y) g"{y) x 



3f(x) 3 - 3f(x)g'(y) 2 - Af(x)f'(x)f"(x) + Ag(y)f'(x)f"(x) + g(yff {3) (x) 



- 2f{x)f{x)g"{y) + 2g{y)f{x)g"{y) + f(x) 2 f'\x) - 2 f(x) g(y) f^\x] 

(10) 

This elimination of the g"'(y) derivatives gives us the equation F(f(x), g(y)) = 
we choose to work with. The factorisation we encounter is one of the simplifica- 
tions we drew attention to earlier. Of course we could have taken the k = 1, / = 2 
and s = 3 form of (||) as our equation F(f(x),g(y)) = 0. It appears that this, 
together with one further partial derivative, d x F(f(x), g{y)) = is sufficient to 
determine g{y), but at the expense of far greater work. In particular there is 
no similar factorisation to that we encountered above. The proof presented was 
devised to circumvent the difficulties of this latter route. 
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Now the first term appearing on the right of (|T0| ) may be written as 

/ f ( \ , \\i d ( fi x ) -9'(y)\ 

(/w - 9(!,)) ^ I J- 

while the second term may be expressed as 

{f{x)-g{y)f d / fjx) (f'{x) 2 -sf{yf) _ 2f'(x) f'jx) f^jx) \ 
9'iv) dy\ (f( x )- g{y )f (f(x)-g(y)) 2 fix) -giy))' 



Our assumption of nonconstant solutions means that fix) ^ g(y) and g\y) ^ 0. 
Thus the vanishing of (|T0|) means either 

!M - 9'iv) _ c (x) 

fix) - giy) ~ Cl[x) 

or 

f'jx) if( X ) 2 -g'(yf) _ 2 f'(x) f{x) f^jx) = 

if(x)-g(y)f (f(x)-g(y)) 2 f{x) - g{y) 2{%h 



according to whether the first or second terms in (|T0D vanish. Therefore, assum- 
ing only third derivatives exist, we have shown (after rearranging) that either 

g'iy) = hg(y) + k, or g\yf = p 3 giy) 3 + p 2 giy) 2 + Pi giy) +p - (11) 



Both are cases of (f|). A consequence of ( |TT| ) is that derivatives to all orders exist 
for solutions of either of these differential equations. The (analytic) solutions 
of these differential equations may be expressed (generically) in terms of the 
exponential and Weierstrass p-function, and indeed the exponential solution of 
the linear differential equation corresponds to a degeneration of the p-function 
differential equation. Further, the identical argument but upon interchanging 
the role of y and x shows that f(x) is also a p-function or a degeneration. We 
have however yet to establish that / and g satisfy the same differential equation. 
Before turning to this there is one point that needs to be clarified. 

In principle it is possible for the function g(y) giving the vanishing of flUf ) to 
be a solution of one of the differential equations (|TTD in one domain and satisfy the 
other differential equation outside of it, with g(y) and its first three derivatives 
matching at any boundary. Such a possibility does not arise in our problem. 
One can readily show that matching a solution g\ (y) of the first differential 
equation ( pT|) with a solution g 2 iy) of the second differential equation at a point 
y , and requiring the first four derivatives to agree at this point, entails gi(y) = 
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Q2{y)- Thus the solutions to the vanishing of (|T0| ) being envisaged perforce have 
discontinuous fourth derivative at such points Uq. Now for our problem, the 
coefficients Zj and Pi of the differential equations (|TTJ) are in fact functions of x, 
given explicitly below. We will shortly see that it is this aspect of our problem 
that rules out the functions g(y) envisaged in this paragraph. 

We now establish that the functions / and g satisfy the same differential 
equation by determining the coefficients U and pi appearing above in two different 
ways. By directly differentiating ([11]) one obtains (for example) that 

, = sTM g'{y)9 {A) {y)-9"{y)g ( " ) {y) 

(The remaining coefficients will be listed below.) Alternately the coefficients 
may be determined by evaluating the functions C\^{x) arising upon integration. 
These functions may be determined in a variety of ways. If (|T(]) vanishes, so 
does the derivative with respect to x of the right-hand side of the equation. 
Substituting either of the expressions for g"{y) corresponding to the two terms 
of fllOD into this derivative equation yield equations for g'(y). Comparison with 
the equations defining C\^{x) show (with a little work) 

rM /"(*) rM f(x)f^(x)-f"(x)f^(x) 



/'W zw 3/'(x) 

(The same expressions arise by applying L'Hospital's rule to the equations defin- 
ing C 1)2 (x).) The coefficients ^ and are then determined to be 



h 



f(xY - f(x) f"(x) 
f'(x) 

f"(x) 
f'(x) 



Po = f{xf - 2f(x) f"(x) + f{x)2 f Q {x) ~ (f(x) /«>(*) - f"(x) f 3 \x)) 

V, = 2 f'{x) - 2fiX ^ {X) + g (/'(x) f^(x) - f"(x) f^(x)) 

P2= iw~iw {nx) /(4)(x) ~ f " {x) /(3)(x)) 

_ f'(x)f^(x)-f"(x)fW(x) 
P3 ~ 3/'(x) 3 

(12) 
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Thus, for example, we have shown that 



, _ g"(y) _ f"(x) 

1 g'(y) /'(*) ' 

and so l\ is in fact a constant. Similarly p 3 is seen above to be a constant. Fur- 
ther, these constants are determined by expressions symmetric in the interchange 
of / and g and the same is true of the remaining coefficients. This symmetry 
means that / and g will satisfy the same differential equation. Also, if we had 
chosen to work with / and h the above argument shows that they satisfy the same 
differential equation. Therefore, each of /, g and h satisfy the same differential 
equation. Finally observe that the coefficient p^, which is constant, involves the 
fourth derivative of the functions satisfying flTT| ) for arbitrary position. This 
expression shows one cannot construct a constant p^ from two (nonconstant) 
solutions of (|TT]) whose first three derivatives agree and whose fourth derivatives 
differ. This rules out having solutions to ( |TT| ) that satisfy the first equation on 
one region and the second equation elsewhere, and having discontinuous fourth 
derivative on the boundaries. 

□ 

Corollary 1 A nonconstant solution f(x) of (|I[) satisfies the differential equa- 
tion (fjp. 

Remarks: (i) If one was happy to assume the five-times differentiability of the 
function f(x) satisfying ([J) then this result may be obtained very quickly. Let 
/ = g in the lemma. We need proceed no further than (^). Simplification arises 
because we have only the one function and its derivatives appearing throughout 
and we may further let y — x yielding (nonlinear) differential equations satisfied 
by f(x). Thus, for example, taking k = 2, I = 4 and upon setting y = x in (|8|) 
we obtain 

o = (fr(Q'(^(Q')'. as) 



r ) \r\r 

(This in fact corresponds to first nonzero term appearing in a Taylor series 
expansion around y = x + e in the k = 1, I = 2 expressions of the lemma.) The 
two final two factors here correspond to the two terms appearing in (|T0|) ; they 
may be straightforwardly integrated to yield the corollary. 
(ii) Another (again not straightforward) route to solving ([!]) is first to express 
this functional equation in a rather different manner. Parameterise x + y + z = 
by x = £ + ? 7) y = £,~V an d z = — 2£. Then with A denoting the standard central 
difference operator (defined by A/(£) = f(C + v)~f(^~ r })) an d H denoting the 
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standard average operator (with = (/(£ + rj) + /(£ — rj))/2) then ([I]) may 

be written 



A/(0 /*/'(£) " A/'(0 nf® = A/(0 F'C*) - A/'(0 F(4 

Here we have set F(z) = /(— 2£) and = — -^F(z)/2 and we may view this 

as an equation in rj with £ fixed. Expanding this equation as a series in 77 leads 
to an infinite set of relations which are first order in F(z), 

a k F\z) -b k F(z) + c k = 0. 

These equations are analogous to those involving a^, foj and q in the lemma. 
Here the coefficients are functions of /(£) and its derivatives. The lowest order 
coefficients are 

ai = /', b 1 = f", Cl = ff"-f 

a 2 = r, b 2 = r, C2 = -4f r + ir + sr 2 . 



Elimination of F(z) amongst these is also sufficient to obtain ([13]). 

At this stage we have found a necessary condition for nonconstant functions 
satisfying (Q) and (|1|): the functions must satisfy (f|) and so are of the form 

ap(5x + 7) + /5, ae* 51 ' + /3, ax + /3. 

As remarked upon in the introduction, the exponential and linear solutions 
clearly satisfy the functional equations without the restriction x + y + z = and 
we need only discuss the first solution here involving the p-function. The noted 
invariances of the functional equations mean we need only determine whether 
any restrictions must be placed on the translation parameter 7 unspecified by 
the differential equation. 

Lemma 2 The function f(x) = p(x + 7) satisfies ([|) provided 37 is a lattice 
point of the p-function. 

Similarly, the functions f(x) = p{x + ^i), g{y) = p{y + ^ 2 ), h(z) = p{z + ^ z ) 
satisfy (^) provided 71 + 72 + 73 is a lattice point of the p-function and these may 
be chosen so that 71 + 72 + 73 = 0. 

Proof: The result follows from the following identity [[F] 

cr(a + b + c) cr(a — b) a(b — c) cr(c — a) 



1 1 1 

p(a) p{b) p(c) 
p\a) p\b) p\c) 



a{afa{bfa{cf 



where a is the Weierstrass sigma function that vanishes at the lattice points 
of the p-function. Letting a = x + 7, b = y + c = z + 7 we see that 
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f(x) = p(x + 7) satisfies (|I|) provided 37 is a lattice point as stated. Similarly 
we see the functions /, g and h given in the lemma satisfy @ provided 71+72+73 
is again a lattice point. Suppose now 71 + 7 2 + 73 = L is such a lattice point, 
and set 7 3 = — (71 + 72)- The periodicity of the jp-function means 

h(z) = p(z + 73) = p(z + 73 - L) = p(z + 73) 

and so we may choose the solutions of @ so that 71 + 72 + 73 = 0. □ 
It remains to discuss the case when at least one of the functions in (^) is a 
constant which we may take to be h(z). 

Lemma 3 Let f, g, h satisfy (j^j. If h(z) is a constant then either 

1. One of the functions f(x) or g(y) is the same constant as h(z), in which 
case the remaining function is arbitrary, or 

2. up to the invariance of f(x) = g(x) = e x . 

Proof: Using the invariance of (0) we may suppose without loss of generality 
that h(z) = 0. Then (for all x, y) 







f(x) g{y) 
f'{x) g'(y) 



(d y -d x )f(x)g(y). (14) 



If neither of f(x) or g(y) vanish identically then f(x) = e Sx+ei and g(y) = e Sy+t2 , 
with ei ( 2 arbitrary and S possibly zero (giving the constant solutions). Using the 
invariance of (|2|) we may choose f{x) = g(x) = e x . Finally if (say) f{x) = we 
see that g(y) is arbitrary. □ 
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